Novel effects of electromagnetic interaction on the correlation of nucleons in nuclear 

matter 
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The electromagnetic(EM) interactions between charged protons on the correlations of nucleons 
are discussed by introducing the Anderson-Higgs mechanism of broken U{1) EM symmetry into 
the relativistic nuclear theory with a parametric photon mass. The non-saturating Coulomb force 
contribution is emphasized on the equation of state of nuclear matter with charge symmetry break- 
ing(CSB) at finite temperature and the breached ^5*0 pairing correlations of proton-proton and 
neutron-neutron. The universal properties given by an order parameter field with a non-zero vacuum 
expectation value (VEV) nearby phase transition are explored within the mean field theory(MFT) 
level. This mechanism can be extended to the charged or charge neutralized strongly coupling 
multi-components system for the discussion of binding or pairing issues. 
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Understanding the properties of nuclear matter un- 
der both normal and extreme conditions is of great im- 
portance in relativistic heavy ion collisions and explain- 
ing the appearance of compact objects such as the neu- 
tron stars and neutron-rich matter or nuclei. The de- 
termination of the properties of nuclear matter as func- 
tions of density/temperature, the ratio of protons to neu- 
trons, and the pairing correlations-superfluidity or super- 
conductivity is a fundamental problem in contemporary 
physics P, The discussion about the property of nu- 
clear ground state-binding energy and pairing correla- 
tions at low temperature is substantial. 

The in-mcdium behavior associated with the many- 
body characteristic is the key, while a non-perturbative 
approach is crucial. The theoretical difficulty of mak- 
ing low energy calculation directly with the fundamental 
quantum chromodynamics(QCD) makes effective theo- 
ries still desirable. As accepted widely, the relativistic 
nuclear theory can successfully describe the saturation at 
normal nuclear density and the spin-orbit splitting 
The further developments oi a — uj theory of quantum 
hadrodynamics model (QHD) make it possible to deter- 
mine the model parameters analytically from a specified 
set of zero-temperature nuclear properties and allow us to 
study the hot nuclear properties and study variations of 
these results to nuclear compressibility or pairing correla- 
tion, and even the symmetry energy coefficient according 
to baryon density which are not well known. Further- 
more, the in- medium hadronic property has attracted 
much attention with this kind of models and there are 
many existed works although the relation between QCD 
and QHD has not been well established^, Q- In recent 
years, with the refinement of nuclear theory study, the 
a — oj theory has been placed in the context of effective 
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theory and it is argued that the vacuum physics has been 
explored in part by this kind of modelS(9..1fl.J. In physics, 
with the obvious non- vanishing fermion nucleon mass in 
relevant Lagrangian, the hidden chiral symmetry is ex- 
plicitly broken. 

The theoretical S-wave pairing correlation issue is a 
long-standing problem. The fundamental ^Sq pairing in 
infinite nuclear matter within the frame of relativistic 
nuclear field theory was first discussed by H. Kucharek 
and P. Ring0|, and it was found that the gaps are 
always larger for three times than the non-relativistic 
results [l^ Il3|. Especially, the very uncomfortable non- 
zero gaps of ^Sq pairing correlation at zero baryon den- 
sity obtained with frozen meson propagators in rela- 
tivistic field theory, as recently pointed out by us[l^. 
remind us that the realistic nuclear ground state with 
MFT approach might not be EM empty. On the other 
hand, the well established low temperature supercon- 
ductivity theory tells us that it would be very inter- 
esting to discuss the broken local EM symmetry effects 
on the properties of the nucleons system. Although the 
in-medium nucleon-nucleon interaction potential induced 
by polarization can give a significant ly im proved descrip- 
tion for EOS and superfluiditv[^ fl5| . the pairing dif- 
ference of PP (proton-proton) from NP (neutron-proton) 
or NN(neutron-neutron) has been discarded. Other ap- 
proaches also recently found that polarization effects sup- 
press the S-wave gaps by a factor of 3 ~ 4^^. The nu- 
merical magnitude of ^5*0 gaps is not sensitive to a special 
parameters set and integral momentum cutoff when the 
polarization effect is taken into account [l^. 

The pure neutron matter cannot exist in nature, and 
the realistic nuclear matter is subject to the long range 
EM interaction. The changes of symmetry properties 
associated with possible phase transition realized on 
some conditions attract physicists very much. In nuclear 
physics, charge symmetry breaking explored by the quite 
different empirical negative scattering lengths ajvAr(p) 
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and app is a fundamental fact[l3 and there are ex- 
isted works to address its theoretical origin^. Coulomb 
correlation effects are a fundamental problem in nuclear 
physics and play an important role for the property of nu- 
clear matter|l9l E3 | , which may lead to rich phase struc- 
tures in the low temperature occasion. For example, in 
Ref. [iJl the influence of the non-saturating Coulomb in- 
teraction is recently incorporated in the multi-canonical 
formalism attempting to explain the reported experimen- 
tal signatures of thermodynamic anomalies and the pos- 
sible liquid-gas(LG) phase transitions of charged atomic 
clusters and nuclei 22] . One may naturally worry about 
the important role of the Coulomb repulsion force on the 
properties of charged system and the thermodynamics 
of charged/neutral nuclear matter to be reflected by rel- 
ativistic nuclear theory and corresponding approaches. 
Within the models based on a — uj field theory and usual 
adopted approaches such as MFT or relativistic Hartree 
approximation(RHA), one can suppose the similar in- 
teractions between PP and NP or NN's, with the weak 
EM interaction being neglected compared to the residual 
strong interaction between nucleons. Theoretically, the 
direct (Hartree) Coulomb contribution of charged pro- 
tons to the EOS cannot be included due to the Furry 
theorem's limit. Although the exchange (Fock) contribu- 
tion can be included in principle from the point of view of 
field theory, the involved calculation and radioactive cor- 
rections caused by relevant infrared singularity of photon 
propagator still remain to be done even in the relatively 
simpler zero-temperature occasion in nuclear physics. If 
one asks what the difference between PP and NP or NN 
pairing correlation is, the original version of QHD with 
MFT or RHA approaches cannot tell us anything. Al- 
though one can expect that the isospin breaking coupling 
terms such as pNN etc. might reflect the Coulomb re- 
pulsion contribution on the thermodynamics of isospin 
asymmetric system to some extent, the pairing differ- 
ences between PP and NN, NP exist even for symmetric 
nuclear matter incorporated with the quite different em- 
pirical negative scattering lengths ajvAr(_p) and app. How 
to incorporate the important role of EM interaction with 
CSB on the thermodynamics of charged/neutral system 
or the property of nuclear ground state on a microscopic 
level (continuum field theory) remains an intriguing task 
even in an oversimplified way (MFT or RHA ) but with 
thermodynamics self-consistency. 

In this Letter, we propose a systematic way to perform 
the link between the bulk and pairing correlation many- 
body properties of charged/neutral two-components nu- 
cleon systems through a relativistic nuclear field theory 
involving the interaction of Dirac nucleons with massive 
photons as well as the well known scalar /vector mesons. 
Inspired by the continuum field theory of phase transition 
and based on QHD-II, the constructed phenomenological 
Proca-like Lagrangian through Anderson-Higgs mecha- 
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where cr, w^, and A^ are the scalar-isoscalar, vector- 
isoscalar, vector-isovector meson fields, EM field with the 
field stresses 

Ffj.!/ = d^An — duApi^ 

for u), p and A^'s, respectively. The M, TOo-, m^, nip and 
are the nucleon, meson and photon masses, while go-, 
<?a;, <?p and e are the coupling constants for corresponding 
Yukawa-like effective interaction, respectively. 

Here the Lagrangian with CSB does not respect the 
local Uem{^) gauge symmetry which is broken by the 
ground state with non-zero local electric charge of pro- 
tons (although the system can be globally neutralized by 
the surrounding such as electrons to maintain the sta- 
bility for compact object through /3-equilibrium) . Also, 
the quartic-cubic terms of a nonlinear-interaction poten- 
tial U{cr) = ba^ + ccr^ with the additional phenomeno- 
logically determined parameters h and c have not been 
obviously preferred in order to discuss in a more general 
way although a specific assumption in U{a) can give a 
reasonable bulk compressibility for nuclear matter. 

The mean field approximation can be used to discuss 
the thermodynamics of charged nuclear matter, from 
which the effective potential is derived in terms of finite 
temperature field theorV|,,1..4. .25.] 

O/TA ^ 2 i2 1 2 2 1 2 2 1 2 /i2 

-T^V / d3k{ln(l + e-^(^*-^»*)) 



(2^) 
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+ ln(f + e-'3(^* +''*))}, (2) 



where i — P, N represents the index of proton (P) and 
neutron (N), respectively and V is the volume of the 
system. With the thermodynamics relation 
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one can obtain the energy density 



Pp 



(3) 



and pressure p = — fl/V. The baryon density is 
PB =< >b=Ypi, 
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In above expressions, ni{fi*,T), ni{p*,T) are the dis- 
tribution functions for (anti-)particles with E* — 
y/\^ + M*^ . The effective nucleon mass M*, chemical 
potentials /^p(jv) introduced by the tadpole diagrams 
of the sigma, omega, rho mesons and photon self-energies, 
respectively 
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where Pp(n) is the proton (neutron) chemical potential. 

The photon mass appears as a free parameter 
which is closely related to the Coulomb energy (reflect- 
ing the binding energy contributed by adding a proton to 
or removing a neutron from the system) and correspond- 
ingly to Coulomb compression modulus Kc- It is worthy 
noting that the Coulomb energy can be discussed by the 
conventional many-body approaches such as the Thomas- 
Fermi theory with variational principle |2fil|. Within rela- 
tivistic MFT, the Kc has been analyzed in the literature 
such as in Ref-ll^l with the scaling model|2^ phenomeno- 
logically. The bulk compression modulus K and Kc are 
defined by 
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Here, pb, Pa, &b £^re the baryon density, the normal 
baryon density and the binding energy per nucleon with 

i?o-[3/(4^Po)]^ K'^ipl- 



'dpi 



The repulsive Coulomb force will modify the EOS sig- 
nificantly for the realistic charged system produced in 
heavy ion collisions. Especially, it can make the criti- 
cal temperature Tc of the LG phase transition decreased 
to a smaller value. With careful numerical study, it 
is found that the softness of bulk EOS (characterized 
by K) is not sensitive to the Coulomb interaction but 
the critical temperature Tc as well as Kc is very sen- 
sitive to this repulsive force. The additional Coulomb 
energy term in pressure and energy density Eq.Q con- 
tributes to removing the theoretical instability in the 
high baryon density region caused by a negative param- 
eters set of b and c in the nonlinear self-interaction term 
U{a) (for obtaining a reasonable compressibility modu- 
lus of bulk EOS). One can estimate the parameter 
is about 20 — 30 MeV for a reasonable critical tempera- 
ture Tc ^ 16 MeV accessible in heavy ion collisions. In 
Fig^ we give the curves of pressure versus baryon den- 
sity and the Coulomb compression modulus Kc accord- 
ing to the interaction strength characterized by with 
frozen parameters go- and g^j which are determined by fit- 
ting the binding energy e° = —15.75 MeV and the bulk 
symmetry energy coefficient asym = 35 MeV (for sym- 
metric nuclear matter at the empirical saturation den- 
sity po = 0.1484/m-3 with T = 0)0]. The qualitative 
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FIG. 1: For charged nuclear matter without considering the 
charge neutral condition with the set (a) of Table|I] (a) Pres- 
sure versus rescaled density ps/po with (solid) and without 
(dashed) the Coulomb repulsion interaction, (b) Kc versus 
the order parameter m-y describing to what extent the EM 
symmetry is broken. 

Coulomb effect on the deformation of phase space distri- 
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bution functions resulting from Eq. Q and Eq. ^ can be 
reflected by the proton fraction ratio: Yp = pp/ps- It 
is found that this ratio changes significantly according to 
temperature T and total baryon density pb- 

Therefore the electric repulsive force plays an isospin 
violating role for the many-body property. Indeed, there 
is some kind competition between the p and photon's 
isospin breaking effect on the phase space distribution 
function deformation. The Coulomb force makes the pro- 
ton fraction decreased while the p meson plays a weak in- 
verse role. Furthermore, one can readily derive the sym- 
metry energy coefficient formula at T = 

asym — 2 Q^2 l*=0 

^2 
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127r^ 
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In fact, if without taking into account the repulsive 
Coulomb contribution, one must introduce a very large 
coupling constant gp to approach the empirical symmet- 
ric coefficient agym which is very far from the empirical 
coupling constant gpNN extracted experimentally. From 
Eq.(jHl, if taking gp = 2.630 and asym = 35 MeV, the 
free parameter can be fixed accordingly. With close 
study, the numerical magnitude of is more sensitive 
to M* (and hence the softness of bulk EOS ) than to gp. 
The relevant parameters are listed in TableClcorrespond- 
ing to the L2 set of Ref.|3] except of taking gp — 2.63. 
Further exact fitting to finite nuclei data such as charge 
density radii distribution etc. 0,113 and addressing the 
spin-orbit splitting issue with the mirror-symmetry topic 
can contribute to giving a solid limit for fixing and 
gp. The large tensor and spin-orbit forces are also crucial 
for understanding finite nuclei and neutron star struc- 
ture which can be explored through the study for mirror- 
nuclei. Let us mention that the hitherto overlooked but 
important EM interaction role on the spin-orbit splitting 
for some mirror-nuclei is recently found in Ref.|30j. 

To study its effect on the electric neutral nuclear mat- 
ter such as compact proto-neutron star would also be 
interesting. For the simplest NP + e + system stabi- 
lized through /3-equilibrium, charge neutrality condition 
makes the proton fraction ratio very small. It is found 
that the Coulomb force does not modify this picture as 
indicated by Fig|21 This is consistent with above result 
that the Coulomb interaction doesn't change the softness 
of bulk EOS significantly. 

Correlations not only do manifest themselves in the 
bulk properties but also modify the quasi-particle prop- 
erties of nucleons in a substantial way. Conceptually, 
the PP and NN(P) pairing correlations should be quite 
different from each other. The former has additional su- 
perconductivity contribution due to the electric charge of 
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FIG. 2: For electric neutral nuclear matter stabilized through 
/3-equilibrium: (a) Pressure versus rescaled density, (b) Pro- 
ton fraction ratio Yp. Line-styles are similar to FiglTI 



protons in addition to the attractive residual strong in- 
teraction compared with the scenario of NN correlation. 
For the fundamental ^5*0 pairing, the energy gap equa- 
tion of nucleon-nucleon pairing in the frame of relativistic 
nuclear theory can be reduced to[ni[ll[l3 
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and the coupled effective mass gap equation has been 
neglected here for brevity. 

In Eq.Q, the assymetrized matrix elements Vpp{p, k) is 
obtained through the integration of v{p, k) over the angle 
9 between the thrcc-momcntums p and k with w(p, k) 
being the particle-particle interaction potential 



«(p,k) = T 



M* 



2E*{k)E*{p) 
Tr[A+(k)rA+(p)7"T+r+T7"] 
(k - p)2 + m% 



where A-|.(k) = ^^ii* projection operator of the 

positive energy solution and T = 17^7'^ is the time re- 
versal operator. The F is the corresponding interaction 
vertex of (J /u>, p{'~i) with nucleons while m-y is the photon 
mass. This static electric contribution to the gap has 
been indicated in FigOl as a curve of gap versus density 
with a integral momentum cutoff = 3.6/to~^ and the 
set (a) in TableQ] to numerically solve the integral gap 
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equation. The difference between PP and NN pairing 
correlations reflects that virtual photons proceeding in 
the space-like momentum transfer regime carry a unique 
information on the EM properties of nucleon interaction 
responsible for the nucleon structure. As indicated by 
Figl^l the long range but screened Coulomb interaction 
affects the correlation function of proton-proton pairing 
significantly esp. in the low momentum regime. 
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FIG. 3: (a) Pairing gap A/ at the Fermi surface versus Fermi 
momentum kf. (b) Gap function A(fc) versus momentum k 
for fixed Fermi momentum. The solid line corresponds to the 
result of proton-proton pairing correlation and dashed line to 
that of neutron-neutron. 

The physical reason for the parametric description 
of the EM interaction in this approach is that at 
first one can note the existence of locally charged sys- 
tem/cluster, i.e., the electromagnetic field condensation 
~< 'ijjp^'^tpp > (corresponding to the spontaneously 
breaking of local gauge symmetry while the gauge field 
obtains mass) in the low energy scale although the stable 
system should be globally neutral with surrounding such 
as electrons through /3-equilibrium. This is very much 
similar to the chiral condensation '--^< ^pip > at low en- 
ergy scale. Second, from the point of view of continuum 
field theory with symmetry changes, the physics back- 
ground of well known low temperature LG phase tran- 
sition still remains to be explored. Especially, how to 
reflect the CSB characteristic in relevant effective theory 
and approaches remains to be performed. Third, in the 
multi-components Fermi/Bose systems the CSB would 
lead to more rich phenomena, e.g., compared with the 
metal electric superconductivity occasion(the ions fixed 
as lattice). 



From the point of view of Maxwell QED, because pho- 
tons are massless, photon-mediated interactions are long 
range in contrast with a point-like meson-nucleon interac- 
tion in the existed QHD-like Lagrangian. The long range 
nature of photon-exchange manifests itself in the infrared 
singular behavior of the photon propagator. This charac- 
teristic enhances the contribution of very soft, coUinear 
photons to the correlation energy for the EOS or the 
pairing problem by noting that this divergence should be 
avoided by the resummation approach as done in QCD or 
QED. Essentially, different from the QCD occasion(with 
the magnetic mass cutoff due to the non-Abelian self- 
interaction of gluonsl [31) . there is no magnetic screen- 
ing in QED, which makes it very involved to discuss the 
superconductive behaviors in strong magnetic field oc- 
casion such as in compact star environment with con- 
ventional QHD-like Lagrangian. This approach makes 
it possible to further study Meissner as well as Debye 
screening effects in such as astrophysics The premise 
of this approach as a nonlocal effective theory nearby a 
phase transition with CSB through one-meson(photon)- 
exchange picture would be very powerful in addressing 
nuclear matter (either symmetric or asymmetric) many- 
body properties and even those of finite nuclei. 

In summary, the long range non-saturating Coulomb 
interaction plays an important role in the property of 
charged/neutral nuclear matter, which can be incorpo- 
rated simultaneously with the residual strong interaction 
within MET of relativistic nuclear theory through an ef- 
fective Proca-like Lagrangian. The deformation of phase- 
space distribution function of nucleons attributed to the 
static electric interaction can manifest itself on the ther- 
modynamics or the property of ground state of charged 
nuclear matter and the quasi-particle spectrum while the 
photon mass rrij controls the strength. Especially, the 
repulsive Coulomb force makes the critical temperature 
lower than the existed theoretical anticipation and con- 
tributes to interpreting the accessible experimental re- 
sults. Furthermore, the breached PP and NN pairing 
correlation strengths open a new window for the study of 
nuclear matter EM property/nuclei structure and would 
lead to rich physical phenomena. From the view of point 
of many-body physics, the low-temperature LG phase 
transition found in heavy ion collisions and the differ- 
ent correlation strengths for PP and NN(P) bound-state 
can be seen as the fingerprint of broken EM symmetry 
within MFT to some extent. Our discussion based on as- 
suming the spontaneously broken EM gauge symmetry 
highlights that the U{1) electric charge symmetry violat- 
ing effects should be taken into account simultaneously 
with the SU{2) isospin breaking effects played by such as 
pNN coupling. The weak coupling interaction is mixed 
with other stronger ones and plays an important role for 
the many-body effects. Especially, the overlooked EM in- 
teraction contribution on the many-body properties such 
as thermodynamics, binding, pairing mechanism etc. of 
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TABLE I: The parameters are with M = 939, rup = 770, 
m„ = 783, m„. = 520 MeV(s) and is in (MeV). Cf = 



Set gl 


35 


glicl) 


m^(Cl) 


M* 
M \P0 


MFT 










a 91.64 


191.05 


6.91 (10.28) 


30.44 (87.27) 


0.540 


b 







28.79 (97.55) 




c 




65.58 (97.55) 


oo or e=0 (0) 




RHA 










d 69.98 


102.76 


6.91 (10.28) 


26.636 (113.96) 


0.731 


e 







25.51 (124.24) 




f 




83.54 (124.24) 


oo or e=0 (0) 





nucleons in nuclear matter should be carefully consid- 
ered from the point of view of continuum field theory. 
Acknowledgments: Ji-sheng Chen acknowledges the 
beneficial discussions with Prof. Hong-an Peng(Peking 
University). This work was supported by NSFC under 
grant Nos. 10175026, 90303007. 



[1] P. Daniellewicz, R. Lacey and W. G. Lynch, Science 298, 
1592 (2002); J.M. Lattimer and M. Prakash, ibid., 304, 
536 (2004). 

[2] D.J. Dean and M. Hjorth- Jensen, Rev. Mod. Phys. 75, 
607 (2003). 

[3] J.D. Walecka, Ann. Phys. (N.Y.) 83, 491 (1974); 
B.D. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 
(1986). 

[4] B.D. Serot and J.D. Walecka, Int. J. Mod. Phys. E 6, 
515 (1997). 

[5] J. Boguta and A.R. Bodmer, Nucl. Phys. A 292, 413 
(1977); J. Zimanyi and S.A. Moszkowski, Phys. Rev. C 
42, 1416 (1990). 

[6] B.-A Li, Phys. Rev. Lett. 85, 4221 (2000). 

[7] For example, H. Shiomi and T. Hatsuda, Phys. Lett. 
B 334, 281 (1994); H.-C. Jean, J. Piekarewicz, and 
A.G. Williams, Phys. Rev. C 49, 1981 (1994). J.-S Chen, 
J.-R Li, and P.-F Zhuang, J. High Energy Phys. 11, 
014 (2002)(nucl-th/0202026); Phys. Rev. C 67, 068202 
(2003). 

[8] J.-S Chen, P.-F Zhuang, and J.-R Li, Phys. Rev. C 68, 
045209 (2003). 

[9] G.E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 
(1991); R. Rapp, R. Machleidt, J.W. Durso, and 
G.E. Brown, ibid., 82, 1827 (1999); G.E. Brown and 
M. Rho, Phys. Rep. 363, 85 (2002); C. Song, ibid., 347, 



289 (2001). 

[10] G. Gelmini and B. Ritzi, Phys. Lett. B 357, 431 (1995); 

G.E. Brown and M. Rho, Nucl. Phys. A 596, 503 (1996). 
[11] H. Kucharek and P. Ring, Z. Phys. A 339, 23 (1991). 
[12] J. Decharge and D. Gogny, Phys. Rev. C 21, 1568 (1980). 
[13] For example, M. Matsuzaki, Phys. Rev. C 58, 3407 

(1998) ; M. Serra, A. Rummel, and P. Ring, ibid., 
65, 014304 (2001); B.V. Carlson, T. Frederico and 

F. B. Guimaraes, ibid., 56, 3097 (1997); B.F. Haas, T. 
Frederico, B.V. Carlson and F.B. Guimaraes, Nucl. Phys. 
A 728, 379 (2003). 

[14] J.-S Chen, P.-F Zhuang, and J.-R Li, Phys. Lett. B 585, 
85 (2004). 

[15] X.-D Ji, Phys. Lett. B 208, 19 (1988); A. Bhattacharyya, 
S. K. Ghosh and S. C. Phatak, Phys. Rev. C 60, 044903 

(1999) . 

[16] For example, A. Schwenk and B. Friman, Phys. Rev. 
Lett. 92, 082501 (2004); A. Schwenk, B. Friman, and 

G. E.Brown, Nucl. Phys. A 713, 191 (2003); ibid., 703, 
745 (2002). 

[17] For example, R.B. Wiringa, V.G.J. Stoks, and R. Schi- 
avilla, Phys. Rev. C 51, 38 (1995); R. Machleidt, ibid., 
63, 024001 (2001); T. Duguet, ibid., 69, 054317 (2004). 

[18] S. Gardner, C.J. Horowitz and J. Piekarewicz, Phys. Rev. 
Lett. 75, 2462 (1995) and references therein. 

[19] A. Bulgac and V.R. Shaginyan, Phys. Lett. B 469, 1 

(1999) ; M. Bender, P.H. Heenen, and P.G. Reinhard, 
Rev. Mod. Phys. 75, 121 (2003). 

[20] T. Maruyama et al, preprint nucl-th/0402002 and refer- 
ences therein. 

[21] F. GulmineUi et al, Phys. Rev. Lett. 91, 202701 (2003). 

[22] P.J. Siemens, Nature (London) 305, 410 (1983); M. 
Schmidt et al, Phys. Rev. Lett. 86, 1191 (2001); F. Go- 
bet et al, ibid., 89, 183403 (2002); S. Siem et al, Phys. 
Rev. C 65, 044318 (2002). 

[23] M.E Peskin and D.V Scroeder, An introduction to Quan- 
tum Field Theory, Weber State University, 1995. 

[24] A. Jenkins, Phys. Rev. D 69, 105007 (2004) and refer- 
ences therein. 

[25] J.I. Kapusta, Finite Temperature Field Theory (Cam- 
bridge University Press, Cambridge, England, 1989) 

[26] For example, K.-S Cheng, C.-C Yao and Z.-G Dai, Phys. 
Rev. C 55, 2092 (1997); H. Shen, ibid., 65, 035802 
(2002). 

[27] H. Kouno et al, Phys. Rev. C 53, 2542 (1996). 
[28] J.-P Blaizot, Phys. Rep. 64, 171 (1980). 
[29] C.J. Horowitz and J. Piekarewicz, Phys. Rev. Lett. 86, 
5647 (2001). 

[30] J. Ekman et al, Phys. Rev. Lett. 92, 132502 (2004). 
[31] For example. Dirk H. Rischke, Phys. Rev. D 62, 054017 

(2000) . 

[32] K.B.W. Buckley, M.A. Metlitski, and A.R. Zhitnitsky 
Phys. Rev. Lett. 92, 151102 (2004); Phys. Rev. C 69, 
055803 (2004). 



